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^ Abstract 

A two-dimensional electron gas exposed to a tilted magnetic field is considered with the Rashba spin-orbit interaction and the 
r*^ ^eeman effect. An exact solution for the eigenvalues was obtained assuming that two opposite spin states of adjacent Landau levels 
D have equal probability. No crossings between adjacent eigenenergies were observed, for the tilt angles studied here {6 < 80°), 

unlike in the perpendicular-magnetic-field case. The absence of crossings lead to quenched beating structures in the oscillations of 
T— I the density of states (DOS). Persistent spin-splittings were observed at the weak magnetic field region. The splittings, however, can 

)}e effectively screened by an increased Landau level broadening. The results shed light on how spins can be controlled through the 

■Rashba interaction strength, the disorder-related broadening and the magnetic field tilt angle. 
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1. Introduction 

The focus and interest in the development of spintronics lies 
in the maneuvering of the spin, rather than solely the charge, 
as the medium of transport in microelectronics lUl IH. For 
instance, the spin field-effect transistor (spin-FET) proposed 
by Datta and Das consists of a two-dimensional electron gas 
(2DEG) connected to ferromagnetic electrodes [2;]. Here spin- 
polarized electrons from a ferromagnetic source electrode are 
injected into the 2DEG. The spin can be rotated because an 
(effective magnetic field, modulated by the gate voltage, is in- 
duced by the spin-orbit interaction. 

Rashba first studied a wurtzite lattice semiconductor where 
an energy extremum is reached not at isolated points in the Bril- 
louin zone but over an extended region |4]. He considered that 
spin-orbit interaction (SOI) serves as a measure of the depth of 
this extremum loop. Using the effective mass approximation, 
Rashba obtained the eigenvalues of electrons with SOI under a 
perpendicular magnetic field. His work has since then served as 
,a benchmark of SOI studies arising from inversion asymmetric 
potentials. 

, Whenever a 2DEG layer is sandwiched in a host system that 
'lacks structure inversion symmetry the spin degeneracy of the 
electron states is lifted even in the absence of an external mag- 
netic field 1^. This zero-magnetic-field spin splitting is rela- 
tivistic in nature. With respect to a laboratory frame, the elec- 
trons are moving in an electric field brought about by the po- 
tential well asymmetry. In the electrons' rest frame, however, it 
is the electric field that is varying. Hence, the electric field ap- 
pears to the electrons as a magnetic field, which is now called 
the Rashba field Br. This nonzero B« has been measured, for 
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example, in GaAs and InSb and other III-V heterostructures and 
in some narrow-gap II- VI semiconductor devices |6, 7, 8, 9]. 

Most research investigations on 2DEG with Rashba SOI have 
dealt with the case when an external magnetic field B is ori- 
ented perpendicular to the plane [IQ, H]. A common feature 
of these studies is the manifestation of beating patterns in the 
Shubnikov-de Haas (SdH) oscillations which is attributed to the 
unequally spaced energy spectrum brought about by the Rashba 
SOI112i]. In general, however, there are two types of spin-orbit 
interactions that govern the nature of a 2DEG found in het- 
erostructures, namely, the Rashba and the Dresselhaus terms. 
The latter's electric field is associated with the host crystal's 
inversion asymmetry and is enhanced by the confinement [5]] 
and by temperature |13]. The Rashba contribution is usually 
stronger as evidenced by III - V semiconductors and also 
in HgTe quantum wells [13]. When only one type of SOI is 
present or when one dominates the other, the beats in the os- 
cillations are expected. But when the two SOI contributions 
are equal, the SdH oscillations are described by a singular fre- 
quency, that is, the beats are suppressed. The absence of the 
beating pattern in the resistivity versus B, for example, is at- 
tributed to equidistant energy levels 

Although most studies dwelt on the perpendicular-magnetic- 
field problem some also explored the effects of 
tilting the orientation of B with respect to the 2DEG plane. In 
an experimental work on an AlGaAs/GaAs heteroj unction, for 
example, it was found that the magnetization versus B- exhibits 
beating patterns at large tilt angles [17]. On the contrary, a col- 
lapse of the ringlike structures in the longitudinal resistivity was 
observed for increasing tilt angles. The ringlike structures orig- 
inate from the crossings of two sets of spin split Landau levels 
from different sub-bands [18]. A theoretical work, on the other 
hand, found that when a 2DEG confined to a parabohc poten- 
tial is under a tilted B, its magnetoresistance exhibit beating 
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patterns lfl9ll . The beats occurred when the cyclotron frequency 
is less than the frequency of the confining potential and ^ is 
nearly parallel to the 2DEG plane. Another research group ob- 
served a beating pattern in the Fermi energy versus the tilt angle 
of a 2DEG under a parabolic potential, free of spin-orbit inter- 
action but inclusive of Zeeman splitting |20]. When the same 
group included the Rashba SOI no pronounced beats were ob- 
served. In between adjacent quantum Hall plateaus, they found 
that in the presence of Rashba SOI or tilted fields, new interme- 
diate plateaus arise. These are attributed to the lifted degener- 
acy ||2ll]. 

The aim of this work is to solve the eigenvalues of a 2DEG 
with Rashba and Zeeman interactions when the external mag- 
netic field is applied at an angle relative to the gas plane's nor- 
mal. No exact solution to this problem is established yet. After 
obtaining the eigenenergies, the density of states will be calcu- 
lated numerically in order to compare the results with experi- 
mental data and other theoretical studies. 



2. Theoretical Formalism: An Analytic Solution 

2.1. The Rashba SOI Hamiltonian 

The Hamiltonian of an electron in the presence of a magnetic 
field B - (Bf, By, B ) can be written as 
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H — Hn + Hr + Hz — 



2m* 



+ a((T X k) ■ z- ft ■ B, 



(1) 



where Hq is the free particle energy, Hr is the Rashba spin- 
orbit interaction, and Hz is the Zeeman energy. Here m* is the 
electron's effective mass. The strength of the Rashba SOI is 
indicated by the parameter a which is assumed to be constant 
for this present work. Here (? are the Pauli matrices and k is 
the wave vector. The magnitude of the latter is determined by 
kj - -Nj + ^Aj, where e is the electronic charge, h is Planck's 
constant h over 2n, and A ^ is the j-th component of the magnetic 
vector potential. 

One expects from Eq. ([TJ that for strong B the Zeeman term 
dominates while the Rashba term leads in the weak field region. 
However, in HgTe quantum wells, SOI is larger or comparable 
to the sp-d exchange-interaction-induced giant Zeeman split- 
ting II13II . Note that the spin-splitting brought about by Hr does 
not require an external B unlike Hz- The magnitude of Hr is 
determined by the potential well asymmetry or by an external 
strain [22] while that of Hz is regulated mainly by the intensity 
of^. 

Electrons in the x - y plane, when subject to a perpendicu- 
lar field of magnitude B, follow cyclotron orbits with angular 
frequency ojc = eB/m*. Hence, they are essentially under a 
harmonic oscillator potential. Here we followed Rashba's for- 
mulation [4] where we let k-t = kj, + iky. When B is tilted at 
an angle 6 with respect to the perpendicular z-axis, the com- 
mutation of the wave vectors depends on the magnitude and 
direction, that is. 



\^kx, ky\ 



ieB, 



ieB 



■ cos 6. 



(2) 



When = 0, we obtain Rashba's results that [k^, A;,,] depends 
on B only. Moreover, since the momentum components do not 
commute they represent incompatible observables with an un- 
derlying uncertainty principle [23]. Relating the wave vector 
to the momentum by — hki, the uncertainty in simultaneous 
measurements of the momentum components becomes 



1 IheB, 



(3) 



Keeping this in mind, the analytic solution can be conve- 
niently obtained by using the ladder operators 



he 



^2eB, 



■k- 



and 



he 



\ 2eB. 



(4) 



which satisfy the canonical commutation relation — 1. 

The Hamiltonian can be rewritten as 



rr _ I Hn + 

' -iTa^+Q.^ 



iTa + D. 
H„-Q., 



(5) 



where H„ - j{aa' + a^'a), Q.j - Bjin* /2B,me is the j-th com- 
ponent of the rescaled Zeeman energy, T - 2a ^f^JhoTc is the 
rescaled Rashba parameter, ^ - m* I2h? and Q.± - Qj- ± /Qy. 
Here we have assumed that the g-factor g - 2. In this work, 
energy units are given in terms of huc where Uc - eB-jm*. We 
used m* - O.OSWf where is the rest mass of an electron. 
In solving the Schrodinger equation. 



H„+Q.^-E if a + Q. 
-iTa' +n. H„-n,-E 



Tjn=0 bnlpn 



= 0, (6) 



we express the wave function solutions in terms of a superpo- 
sition of the harmonic oscillator basis functions (p„ where the 
coefficients a„ and b„ are the unknown spin-up and spin-down 
complex coefficients of the nth Landau level, respectively. This 
yields a set of secular equations 



H- i + - fij a„ + iT ^In + 1 bn+i + 0-b„ = 0, 



(7) 



- ifyfUan-i + n+a„ + \n + - - D., - Ejb„ = 0. (8) 

Note that the secular equations derived by Rashba are re- 
covered when S is normal to the 2DEG plane, that is, when we 
let = f2_ = 0. We will impose this latter condition on sec- 
ular Eqs. O and (O to obtain the 2DEG behavior when 6-0. 
Later we will be needing such information for comparison with 
our results on tilted fields. 

When we consider the pure Zeeman case ( Y - 0), the eigen- 
values are shown to be 



1 

£7 = n + - ± fi. 

^ 2 



(9) 



Accounting for the contribution of the in-plane magnetic field, 

the general Zeeman term Q. - ^Q^. -i- f2j -i- Q? replaces the 

in the perpendicular-magnetic-field case. As expected, the spin- 
splitting brought about by the Zeeman effect Az = 2Q is deter- 
mined by the magnitude of B which includes both the in-plane 
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and perpendicular components. On the other hand, when only 
the Rashba SOI is the source of spin splitting (Q. = 0), the 
eigenvalues can be expressed as 



= n + - VTT 



4T2n. 



(10) 



The spin splitting brought about solely by the Rashba SOI 



Ar = 



E'l^ - Ej^ - Vl + 4T^n is directly proportional to the 



strength of the spin-orbit interaction and to the square root of 
Landau level index. The number of filled Landau levels n is 
inversely proportional to = Bcosd. This means that for a 
given direction 6, the Rashba splitting is more pronounced as 
B — > 0. When the intensity B is fixed, Ar becomes larger as 6 
approaches 90°. 

2.2. The Energy Spectrum 

We need to solve Eq. (|7]l and Eq. (H) in order to obtain the 
general eigenvalues when both the Zeeman and Rashba terms 
are present. Here one needs to have additional specifications 
about the 2DEG since there are more unknown coefficients than 
there are secular equations. This makes the tilted-field case dif- 
ficult to deal with. A probable additional equation would be 
the normalization condition but establishing it is not straight- 
forward because of the coupling of adjacent orthogonal states. 
Due to these chal leng es, some researchers resorted to using per- 
turbation theory [24], the continued fraction numerical method 
1 2l|] , and a quasi-classical approach |25.] . But as far as we know, 
an exact solution for the tilted-magnetic-field case has not yet 
been done. In this work we will take advantage of Landau level 
crossings in order to obtain an analytical solution to the prob- 
lem. 

When B is oriented normal to the 2DEG plane, energy level 
crossings are observed. These crossings are caused by the 
widening energy gap A between opposite spins of the same Lan- 
dau level n. The gap widens due either to an increasing Rashba 
SOI strength lllSi |24] or magnetic field 1.261 until opposite spins 
of different Landau levels intersect, that is, E„-ii crosses with 
E„i while E„+ii crosses with E„i. The state \n - 1, T) becomes 
closer to \n, i) than to |n - 1, X). We invoke that \n - 1, t) and 
|n, i) are equally probable. This also applies to \n + and 
|«, t)- Since the coefficients a„ and b„ measure the probabil- 
ity of obtaining the state <p„, we hypothesize that a particular 
solution can be based on the following 
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Using Eq. (fTST i. the secular Eq. (|7]i becomes 

^n + ^+n,- E^ a„ + T Vn + la„ + Q.-bn = 0. 

Similarly, substituting Eq. ( fTTT i. Eq. (jSj becomes 
T ^Jnbn + Q+a„ + [n +]^- Q., - E^b„ ^ Q. 



(11) 
(12) 
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Figure 1 : The eigenvalues as a function of the Rashba pai'ameter. Here B = 0. 1 
T. The main frame corresponds to 6 = 45° while the inset to 9 = 0. 



Therefore, we obtain 



(15) 



Equation ( fTSl l is valid in the tilted-magnetic-field case only. 
Otherwise we need to set Q+ = Q_ = in Eqs. dTji and (ID. 
Substituting Eq. ( fTSl l into Eq. ( fT3] l, we have 

- £ + T V« + l) - £ + T Vn) - = 0, (16) 



where {E^)z 



'^'^ - « -I- 1 ± Q, which is the z-component of Eq. (|9]l- 



This system of equations leads to a quadratic equation in E 
which can be expressed as 

£- - |2 |« -H i j + Y ( Vn + Vn + l)| £ + ^« + ^ j 

+ T2VnVn + 1 + ^« + ijT(Vn+ Vn + l) 

- Q? + aT ( Vn - Vn + l) - 0+Q_ = 0, (17) 
with the following roots 

E - n + - + — i^yfn + V«TT) 

+ i -^[t ( V^T+T - Vn) + 2Q.]^ -I- 4Q+0_. (18) 

The effect of tilting, based on Eq. ( fTSl l. is to increase the gap 
between the spin up and spin down states. This results into 
the suppression of the crossings (See Fig.[T]) that are observed 
when B is perpendicular to the 2DEG plane (See inset.). This 
is true for the range of the Rashba strength studied in Fig. [T] 
and for the case when Q < 9 < 80°. In the presence of an in- 
plane component of the magnetic field B\\, the Landau levels 
are now substantially separated as shown in the main frame of 
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B(T) 

Figure 2: The density of states as a function of B for E x 13.18 meV and for 
two different Rashba parameters. Note that the two curves coincide at B > 4.5 
T. The plot is taken for T = 0.5 meV and 6 = 22.5°. 

the same figure. The same absence of crossings due to B\\ was 
also shown in Ref. |25] where they considered a quasi-classical 
limit. However, our preliminary investigation shows that be- 
yond 80° crossings reappear. 

Inspecting Fig. [T] closely, for a given n the splitting A is min- 
imal, but not zero, when a = 0. At this point the splitting origi- 
nates solely from the Zeeman term and the energies of the |n, t) 
and \n, i) are close to each other As the Rashba parameter is 
increased, the gap widens until we can observe that the gap be- 
tween |« -I- 1 , X) and |n, T) becomes smaller than that of the |«, t) 
and \n, J,). It is remarkable that the \n+ \,\,) and |«, t) indicate 
parallelism assuring that crossings are unlikely to occur even at 
strong Rashba SOI. 

3. Density of States 
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Figure 3: The density of states as a function of B for E x 13.18 meV in the 
absence and presence of Rashba SOI. Here 8 = 80° and F = 0. 1 meV. 
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Figure 4: The density of states as a function of B for E x 8.79 meV for different 
Landau level widths. Here a = 5.5 X 10"" eVm and 9 = 22.5°. 



In analyzing the properties of a 2DEG we usually refer to 
the density of states (DOS) which indicate the number of one- 
electron levels in the energy range E to E + dE IztIi . In this 
work, we will be employing its Gaussian form 



DOS(£) 



eBcosO 
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{E-Enf 


[in, 




2r2 



(19) 



The tilt angle 9 affects both the magnitude and the center E„ of 
the Gaussian peaks. The parameter F denotes the Landau level 
broadening width. Any variations of the DOS with reference to 
the different tunable parameters involved can be used to predict 
and explain the behavior of thermodynamic quantities. 

3.1. Rashba interaction effects 

For a particular energy level, we plot the DOS versus B for 
two different Rashba SOI strengths in Figs.|2]and|3] In addition 
to the Zeeman term, the presence of the Rashba interaction usu- 
ally widens the spin-splitting gap as can be seen in the eigenval- 
ues obtained in Eq. (fTST l. This is evident in Fig.|2]especially in 



the region 2 T < Z? < 4 T. We can attribute the splitting for dif- 
ferent a to the Rashba interaction exclusively since both curves 
are under the same magnetic field region. At these moderately 
strong fields, the Zeeman splitting is hardly noticeable. In con- 
trast, at a very strong magnetic field region {B > 10 T), changes 
in the Rashba interaction strength will have a negligible effect 
on the DOS. This is expected since the Zeeman term is dom- 
inant at large B. The precise coincidence of the two curves at 
very strong B region supports this. This is further emphasized 
in Fig. [3] where the Rashba term was completely turned off for 
the solid line. At the main frame, the splitting at B > 2 T orig- 
inates from the Zeeman term. From the inset, we notice that at 
weak magnetic fields (B < 1 T, see inset of Fig. O, the spin- 
splitting brought about by the Rashba SOI persists while that 
arising from the Zeeman interaction becomes negligible. 

3.2. Landau level broadening effects 

The broadening parameter is considered to be a measure of 
disorder in the system 112 811 . Based on Fig. |4] wider broaden- 
ing obscures the oscillations and spin splittings at the lower B 
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Figure 5: The density of states as a function of B for £ w 13.18 meV at different 
tilt angles when the Zeeman term is neglected. Here a = 5.5 X 10"" eVm and 
r = 0.1meV. 




Figure 6: The density of states as a function of 6 for E k 8.79 meV showing 
the comparison of the no tilt scenario from the tilted case. Here a = 5.5 X 10"" 
eVm and F = 0. 1 meV. 



region. Together with Fig. |2] we can infer that both the re- 
duction of disorder in the system (smaller F) and increase of 
the heterostructure asymmetry (larger a) pave the way for the 
oscillations at the weak field region to manifest. A similar ob- 
scuring effect was observed for the spinless case ll29ll . Hence, 
samples with smaller F will likely be able to exhibit the pro- 
nounced splittings. 



5.5. Tilt effects 

The presence of an in-plane component of B causes a non- 
uniform phase shift to the DOS oscillations as can be seen in 
Fig. |5] From Eq. ( fT9] l. the height of the DOS peaks is deter- 
mined only by the prefactor B, - BcosO. This explains the taller 
peaks for 6 = 22.5°. The "phase shift" in the DOS oscillations 
that depends on the direction of the field might be related to 
the anisotropy measured in the magnetization data of Ref . iITtIi . 
Since the magnetization leads to a measure of the ground state 



DOS, it is likely that the behavior of the experimental results 
for magnetization reflects a similar shifting of the DOS phase. 

In order to compare the influence of the tilt angle, we study 
Fig.|6] For the perpendicular case (0 = 0), the presence of taller 
peaks in between shorter ones is reminiscent of beating patterns 
observed in thermodynamic quantities Isl |3T|] . The mono- 
tonic increase in the height of the DOS peaks for 6 = 45° can 
point to obscured beats. We can deduce here that the occurence 
of beats can be traced from the crossings of energy levels which 
are present only in the perpendicular-magnetic-field case. This 
is unlike the data of Ref. llTIl where the beat structures become 
more prominent at large tilt angles. We note, however, that the 
measurements reported by their group was with respect to B~ 
only and not with respect to the total B as is the case in this 
present work. Moreover, they have fitted their empirical data 
with a broadening parameter F oc ^/B while a constant F is 
assumed herein. In comparison we also simulated the DOS us- 
ing the VS-dependent broadening and found stronger oscilla- 
tions when B < 1 T. Smeared and broader Gaussian peaks are 
observed when B > 1 T. This is expected because narrower 
widths would be anticipated in weak fields and wider widths 
in stronger fields. Nonetheless there were still no noticeable 
beating patterns. 

4. Conclusions 

A two-dimensional electron gas under a tilted magnetic field 
with Rashba and Zeeman interactions was studied. An exact 
solution was obtained for the case when opposite spin states 
of neighboring Landau levels are equally probable. Eigen- 
value crossings were not observed for the tilt angles studied 
here (6 < 80°). This agrees with the numerical result of 
Ref. [26] wherein the energy crossing degeneracies are lifted 
for a 2DEG without disorder The tilting of magnetic fields for 
(0 < 9 < 80°) resulted to quenching of the beating patterns 
in the density of states. We conclude that it is not the unequal 
spacing between energy levels but the presence of crossings that 
directly cause the beats observed in Uterature. The persistence 
of the spin-splitting in weak and moderately strong magnetic 
fields brought about by the Rashba interaction is more percep- 
tible in systems with less disorder Although the Rashba in- 
teraction causes splittings, it also reduces the amplitude of the 
oscillations similar to the effect of the Landau level broadening. 
The Zeeman contribution to the splitting is comparable to that 
of the Rashba term at moderate fields and becomes dominant at 
strong fields. The determination of these different regimes as 
well as the effects of the field's tilt angle will enable the fitting 
of the 2DEG according to their desired state and application. 

This work shows that tuning the Rashba interaction with 
tilted fields yields 2DEG behavior which varies differently if the 
magnetic field is purely perpendicular. While this manuscript 
was in preparation, the authors came across new results in 
Ref. 13211 which reveal a drastic change in the behavior of the 
longitudinal resistance from the perpendicular to the case with 
tilted fields, even for small angles. How this transition is ap- 
proached must come from self-organizing electron interactions 
in addition to the Rashba and Zeeman interactions. 
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